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Abstract.  Simulations  are  performed  using  a  multiple  temperature  gas  model  to  investigate  translational 
non-equilibrium  effects  in  a  rapid  expansion  of  a  high  temperature  argon  gas  cloud  into  a  rarefied 
atmosphere.  A  set  of  continuum  conservation  equations  based  on  kinetic  theory,  which  includes  anisotropic 
forms  for  the  temperature,  pressure  and  speed  ratio,  are  solved  numerically  using  a  flux  splitting  scheme  for 
the  inviscid  fluxes  and  a  central  difference  scheme  for  the  viscous  fluxes  in  a  time  accurate  manner.  Results 
obtained  for  the  initial  expansion  of  the  spherical  gas  cloud  from  a  high  density  source  condition  show  that 
translational  non-equilibrium  exists  in  the  shock  front  region  which  propagates  into  the  ambient  atmosphere. 
For  a  lower  density  source  condition,  translational  non-equilibrium  not  only  exists  in  the  shock  front  but  also 
in  the  inner  gas  cloud  region  where  the  temperature  normal  to  the  radial  direction  freezes  at  a  value  just  below 
the  initial  source  temperature. 


Introduction 

The  expansion  of  a  high  temperature  gas  cloud  in  a  rarefied  environment  is  a  phenomenon 
that  has  many  applications  in  space  from  supernova  explosions  to  recent  active  space  experiments' 
where  the  interaction  of  an  artificial  plasma  gas  cloud  with  the  ambient  atmosphere  was  studied. 
In  the  rapid  expansion  of  a  gas  cloud,  a  spherical  shock  wave  front  propagates  into  the  ambient 
atmosphere  followed  by  a  contact  surface,  which  separates  the  expanding  gas  cloud  mass  from 
that  of  the  ambient  environment,  and  a  rarefaction  wave  propagates  towards  the  center  of  the 
source  in  the  expanding  gas.  Recent  time-accurate  continuum  simulations^’^  of  a  high  temperature 
metallic  gas  cloud  at  high  altitude  with  radiation  transport  show  that  extremely  high  velocities 
along  with  material  phase  change  can  occur  in  the  flow  of  such  a  gas  cloud  expansion.  Calculated 
values^’^  of  the  continuum  break  down  parameter  proposed  by  Bird  suggest  that  in  some  parts  of 
the  flow,  especially  in  the  shock  region,  the  non-equilibrium  effects  dominate  the  flow  structure. 
Thus,  it  appears  that  particle  simulation  methods  such  as  the  direct  simulation  Monte  Carlo 
(DSMC)  method**  would  be  more  suitable  for  simulating  these  flows.  However,  the  presence  of 
radiation  and  unsteady  behavior  in  these  flows  makes  the  particle  simulations  a  formidable 
computational  task.  Moreover,  the  high  densities  in  the  source  region  along  with  the  phase 
changes  encountered  in  the  flow  present  further  problems  with  the  dilute  gas  assumptions  made  in 
these  particle  simulation  methods.  Experimental  investigation  of  the  initial  evolution  of  such  a 
gas  cloud  is  also  complicated  by  the  very  small  time  scales  involved  and  optical  opaqueness  near 
the  gas  source.  This  suggests  that,  at  present,  a  continuum  method  which  to  some  extent  accounts 
for  non-equilibrium  effects  in  the  flow  will  be  optimum  in  predicting  the  early  expansion  of  such 
a  high  temperature  gas  cloud. 

In  the  present  study  the  time  accurate  numerical  simulations  of  a  rapidly  expanding  high 
temperature  argon  gas  cloud  are  performed  using  continuum  conservation  equations  with 
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anisotropic  temperature  gas  models^"’  This  model  accounts  for  translational  non-equilibrium 
effects  and  is  appropriate  for  these  simulations.  The  equations  are  solved  numerically  using  the 
kinetic  flux  splitting  method*  for  the  inviscid  fluxes  and  central  difference  scheme  for  the 
dissipative  fluxes.  The  purpose  of  this  study  is  to  provide  fundamental  understanding  of  the  initial 
evolution  of  a  rapidly  expanding  gas  cloud  and  propagation  of  the  associated  shock  front  in  the 
ambient  atmosphere  where  translational  non-equilibrium  effects  become  important.  Results  using 
Navier-Stokes  calculations  for  the  dense  source  condition  are  presented  for  comparison. 

Gas  Dynamics  Model 

The  flow  is  governed  by  one-dimensional/spherically  symmetric  conservation  equations  of 
gas  dynamics,  modified  to  account  for  translational  non-equilibrium.  These  modified  equations 
are  derived  from  a  first-order  Chapman-Enskog  expansion  of  the  Boltzmann  equation  using  an 
anisotropic  distribution  function.  This  approach  is  quite  similar  to  that  followed  by  Candler, 
Nijhawan  et  al.m  their  investigation  of  the  normal  shock  structure’.  In  integral- vector  form,  the 
equations  can  be  written  as 
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where  the  conserved  quantities  per  unit  volume  are  given  by 

'P 

P  V 

Q=  "  ■ 

P^tr 

P  ^tp 


Only  one  of  the  two  components  of  the  translational  specific  energy  that  is  perpendicular  to 
the  radial  direction  (e,p)  is  considered  because  the  two  perpendicular  components  are  equal.  The 
gas  dynamics  flux  vectors  are  given  by 
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The  source  terms  considered  in  the  present  formulation  are  assumed  to  be  for  Maxwell 
molecules  gas.  Thus,  fhe  source  terms  are  as  (Ref.7) 
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where  p  is  the  mass  density,  Vr  is  the  radial  velocity,  Tr  is  the  gas  temperature  in  radial  direction, 
Tp  is  the  temperature  in  one  perpendicular  direction  and  R  is  the  gas  constant.  The  expression  for 
krr  and  kpr  above  are  for  Maxwell  molecules. 

Computational  Approach 

The  gas  dynamic  equations  listed  above  are  solved  numerically  in  a  time  accurate  manner. 
The  inviscid  part  of  the  flux,  Fq,  is  calculated  at  each  cell  interface  using  a  kinetic  flux  splitting 
method.*  A  local  ellipsoidal  distribution  function  of  the  form 
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is  then  defined  for  each  cell  and  used  to  formulate  analytical  expressions  for  the  inviscid  fluxes  in 
both  the  radial  and  perpendicular  directions.  The  c  V,  c’p,  and  c’,  are  the  components  of  the  random 
molecular  velocity  in  the  above  expressions.  As  with  other  flux- splitting  techniques,  this  inviscid 
flux  is  formally  first-order  accurate  in  space.  Second-order  accuracy  in  space  is  achieved  by  using 
a  MUSCL  variable-extrapolation  scheme®  in  conjunction  with  a  minmod  slope  limiter  to  ensure 
monotonicity  Viscous  fluxes  are  calculated  with  a  second-order  accurate  central  difference 
scheme. 

Fourth-order  accuracy  in  time  is  obtained  using  a  Runge-Kutta  time  integration  method^®. 
The  time  step  is  chosen  as  the  minimum  of  a  Courant -based  time  step  and  a  step  based  on  the 
magnitude  of  energy  exchange  between  radial  and  perpendicular  directions. 


Results  and  Discussion 


The  high  temperature  argon  gas  cloud  assuming  the  interatomic  potential  to  be  Maxellian 
(inverse  fourth  power)  =  213  K  and  =  2.5 17x10  ^  kg  I  m  s)  is  initialized  as  a  sphere  of 

radius  0.03  m  with  a  uniform  density  and  temperature.  The  initial  density  and  temperature  of  the 
gas  cloud  are  assumed  to  be  20  kg/m^  and  10,000  K,  respectively.  The  argon  cloud  in  equilibrium 
at  this  temperature  would  be  partially  ionized,  but  its  thermal  energy  would  still  be  much  larger 
than  the  electromagnetic  potential  energy.  Thus,  the  electromagnetic  effects  can  be  neglected  for 
the  initial  evolution  of  the  gas  cloud  expansion. 

The  initial  spherical  volume  is  divided  into  240  uniformly  distributed  cells.  One  cell  is 
added  to  the  outer  boundary  after  each  time  step  allowing  grid  to  grow  in  time  as  the  gas  cloud 
expands.  The  procedure  reduces  the  computational  task  during  early  phases  of  the  solution 
because  computation  may  need  up  to  100,000  cells  at  later  times.  The  outer  boundary  condition  is 
assumed  to  be  the  undisturbed  ambient  conditions.  The  ambient  conditions  for  this  study 
correspond  to  60  km  altitude.  Thus,  the  density  and  temperature  used  are  3.09  x  10'“*  kg/m^  and 
247  K,  respectively.  To  focus  on  the  cloud  expansion  and  not  on  molecular  diffusion  and  other 
multi-species  issues,  the  ambient  atmosphere  is  simulated  as  pure  argon  at  the  ambient 
atmospheric  density  and  temperature. 

Figures  1-4  show  the  flow  structure  of  the  spherical  gas  cloud  expansion  at  three  different 
times.  These  figures  show  that  a  well-defined  shock  wave  propagates  outward  in  the  ambient 
atmosphere.  The  shock  is  followed  by  a  contact  surface,  which  separates  the  expanding  gas  from 
the  ambient  environment.  A  strong  rarefaction  wave  travels  behind  the  contact  surface  towards  the 
center  of  the  source.  Figures  1-3  also  show  that  density,  temperature  and  pressure  decrease  in  the 
source  region  as  time  increases.  The  magnitude  of  the  density  ratio  across  the  shock  remains 
almost  the  same  as  the  gas  expands  (Fig.  1).  However,  the  temperature  and  pressure  ratios  across 
the  shock  decrease  gradually  as  the  gas  cloud  expands  (Figs.  2-3). 

Figure  2  shows  that  the  minimum  radial  and  perpendicular  temperatures  in  the  gas  cloud 
behind  shock  front  can  fall  significantly  lower  than  the  ambient  temperature.  This  is  due  to  the 
strong  rarefaction  wave  that  propagates  in  the  opposite  direction  of  the  shock  wave.  Similarly,  the 
radial  and  perpendicular  components  of  the  pressure  behind  the  shock  front  fall  lower  than  the 
ambient  pressure  value  with  time  as  gas  expands  (Fig.  3). 

In  Fig.  4  the  radial  velocity  profiles  are  shown  as  a  function  of  normalized  radial  distance  at 
different  times.  The  radial  velocity  increases  linearly  through  the  cloud  with  increasing  radius  for 
most  of  the  flowfield.  Also,  the  radial  velocity  gradient  decreases  with  time  as  the  gas  cloud 
expands.  The  peak  velocities  predicted  for  the  first  1.2  x  10'“*  s  in  the  gas  cloud  expansion  are  more 
than  3  km/s  (Fig.  4). 

Translational  non-equilibrium  is  evident  in  the  shock  wave  as  the  gas  cloud  expands  (Figs. 
5-6).  Figure  2  and  Figs.  5-6  show  that  the  most  of  the  inner  part  of  the  gas  cloud  remains  in 
translational  equilibrium.  This  results  because  the  mass  density  is  still  high  enough  in  the  cloud  to 
equilibrate  the  translational  modes  at  this  initial  phase  of  the  gas  expansion.  There  is  an  observed 
overshoot  in  the  radial  temperature  just  behind  the  shock  wave,  but  no  overshoot  is  observed  in 
the  perpendicular  temperature  at  these  conditions  (Fig.2). 


Fig  1:  Density  Profiles 


Fig  2:  Temperature  Profiles 


Fig  3:  Pressure  Profiles 


Fig  4:  Radial  Velocity  Profiles 


Figures  5-8  show  the  detailed  shock  wave  front  stracture.  The  radial  and  perpendicular 
temperatures  profiles  near  the  shock  wave  for  two  different  times  (4  x  10'^  s  and  12  x  10'^  s)  are 
shown  in  Figs.  5-6.  It  can  be  seen  from  these  figures  that  thermal  non-equilibrium  exists  just 
behind  the  shock  where  the  radial  and  perpendicular  temperature  differ.  The  radial  temperature 
overshoot  increases  with  time  as  is  demonstrated  in  Figs.  56  where  the  overshoot  in  radial 
temperature  at  2  x  10'^  s  is  about  12%  whereas  at  12  x  10'^  s  the  overshoot  is  about  19%. 
Although  the  ratio  of  the  overshoot  to  the  equilbrium  value  increases  with  time,  the  magnitude 
of  the  overshoot  radial  temperature  decreases  with  time.  The  translational  non-equilibrium  in  the 
shock  wave  region  can  also  be  seen  in  pressure  profiles  (Figs.  7-8),  where,  the  perpendicular 
pressure  in  the  shock  wave  region  lags  behind  the  radial  pressure.  There  is  no  overshoot  in  radial 
or  perpendicular  pressure  profiles  like  those  observed  in  the  radial  temperature  profiles. 


Fig  5:  Shock  Front  Temperature  Profiles  Fig  6:  Shock  Front  Temperature  Profiles 


Fig  7:  Shock  Front  Pressure  Profiles 


Fig  8:  Shock  Front  Pressure  Profiles 


Figure  9  shows  that  the  shock  front  location  from  the  Navier-Stokes  calculations  leads  the 
shock  front  location  from  the  multi-temperature  model  calculations.  This  is  probably  due  to 
differences  in  the  calculation  of  the  thermal  radial  momentum  flux.  For  the  multi-temperature 
model,  the  thermal  radial  momentum  flux  is  equal  to  radial  pressure,  whereas,  for  the  Navier- 
Stokes  equations  it  is  the  sum  of  the  isotropic  pressure  and  the  viscous  stress  in  the  radial 
direction,  which  depends  on  the  radial  velocity  gradient.  In  the  dense  region  of  the  source  the 
density  profiles  from  both  methods  compare  well  (Fig.  9). 

A  second  set  of  calculations  is  also  performed  for  a  rarefied  source  condition.  The  radius  of 
the  initial  source  considered  for  this  set  of  calculations  is  0.05  m  and  the  source  density  and 
temperature  are  100  times  higher  than  the  ambient  values.  The  ambient  values  of  density  and 
temperature  are  7.971  x  10'^  kg/m^  and  194.3  K,  respectively.  These  correspond  to  an  altitude  of 
approximately  100  km.  Figure  10  shows  the  normalized  plots  of  density,  radial  temperature,  and 
normal  temperature  as  a  function  of  the  normalized  radial  distance.  Strong  translational  non¬ 
equilibrium  is  observed  well  into  the  source  region  where  the  perpendicular  temperature  is  frozen 
at  a  value  slightly  lower  than  the  source  value,  whereas  the  radial  temperature  decreases  rapidly 


with  the  expansion  of  the  gas  cloud.  The  perpendicular  temperature  lags  behind  the  rise  in  the 
radial  temperature  in  the  shock  front.  However,  in  the  source  region,  the  radial  temperature  is 
lower  than  the  perpendicular  temperature  due  to  the  large  rarefaction  effects  at  these  conditions. 
Thus,  significant  thermal  non-equilibrium  can  also  be  observed  in  the  source  region  in  the  initial 
phase  of  the  expansion  if  the  initial  density  of  the  source  is  sufficiently  low. 


Fig  9:  Density  Profiles  Comparison  Fig  10:  Density  and  Temperature  Profiles 

Conclusions 

A  kinetic  theory  based  flux  vector  splitting  scheme  has  been  successfully  employed  for  the 
calculations  of  the  inviscid  fluxes  of  the  continuum  model  with  translational  non-equilibrium. 
This  approach  allows  the  analytical  formulation  of  the  inviscid  fluxes  of  the  conserved  quantities 
in  the  different  characteristic  directions. 

Computations  carried  out  for  two  distinct  source  conditions  indicate  that  the  shock  wave 
front  that  propagates  through  an  ambient  atmosphere  will  exhibit  translational  non-equilibrium 
for  all  the  source  conditions  considered.  For  the  less  dense  source  condition,  the  translational 
non-equilibrium  not  only  exists  in  the  region  of  the  shock  front  but  throughout  most  of  the  gas 
cloud  expansion.  The  multi-temperature  model  predicts  a  radial  temperature  overshoot  in  the 
shock  front  for  the  dense  source  conditions  and  this  radial  temperature  overshoot  to  the 
equilibrium  temperature  in  the  shock  increases  as  the  gas  cloud  expands.  Comparison  of  the 
multi-temperature  model  results  with  the  Navier-Stokes  calculations  shows  that  the  shock  front 
location  in  the  Navier-Stokes  calculations  leads  that  predicted  by  the  multi-temperature  model 
because  the  Navier-Stokes  equations  do  not  have  the  translational  non-equilibrium  included  in  the 
formulation  of  the  momentum  flux.  However,  in  the  other  regions  of  the  gas  cloud  expansion 
where  flow  exhibits  translational  equilibrium,  the  two  methods  give  similar  results. 
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